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Abstract 

Take N sites distributed randomly and uniformly on a smooth closed surface. We 
express the expected distance (Dk(N)) from an arbitrary point on the surface to its 
fcth-nearest neighboring site, in terms of the function A(l) giving the area of a disc of 
radius I about that point. We then find two universalities. First, for a flat surface, 
where A(l) = nl 2 , the /c-dependence and the iV-dependence separate in (Dk(N)). All 
/cth-nearest neighbor distances thus have the same scaling law in N. Second, for a 
curved surface, the average f(D^(N)) dfi over the surface is a topological invariant at 
leading and subleading order in a large ./V expansion. The l/N scaling series then 
depends, up through 0(1/N), only on the surface's topology and not on its precise 
shape. We discuss the case of higher dimensions (d > 2), and also interpret our results 
using Regge calculus. 
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1 Introduction 



Many problems arising in applied mathematics involve the distance between neighboring 
sites in a space. One frequently wishes to calculate distances to a nearest neighbor, second- 
nearest neighbor, and more generally, fcth-nearest neighbor. Examples in computational 
geometry and optimization abound, ranging from random packing of spheres to minimum 
spanning trees. Such problems also occur naturally in physics, and have been considered in 
applications ranging from stellar dynamics [|TJ, to interactions in liquid systems, to cellular 
objects such as foams and random lattices 0. 

Here we consider the case of N sites placed randomly, with a uniform distribution, on a 
2-D surface of fixed area. Let the random variable Dk(N) represent the distance between 
a given point x and its /cth-nearest site. The expectation value (Dk(N)) taken over the 
ensemble of randomly placed sites — and in fact all moments (D%(N)) — then exhibit some 
surprising properties. First of all, when the surface is flat (Dk(N)) may be written, up to 
corrections exponentially small in N, as 

. , . . 1 (A; — 1/2)! N\ 

This means that the fc-dependence and the iV-dependence separate; the large N scaling law 
for (Dk{N)) is independent of k. Geometrically, the meaning of this universality is far from 
obvious. Furthermore, the property is not restricted to two dimensions, and turns out to be 
equally valid for flat spaces of any dimension. Second of all, when the surface is curved, while 
this ^-independence no longer holds, one finds on the other hand another universality: if 
(Dk(N)) is written in terms of a l/N series, and averaged over the entire surface, the leading 
and sub leading coefficients of the expansion give topological invariants. Thus to 0(1/N), 
the large N scaling law depends not on the detailed shape of the surface but only on the 
surface's genus. 

In this paper we explore these universalities. We start by expressing (Dk(N)) in terms 
of the area, A(l), of a disc of radius / on an arbitrary surface. We observe that in the 
special case where A{1) consists only of a power of /, the scaling law for (Dk{N)) exhibits 
the universality in k, i.e., is independent of k. Then, we give a relation between A(l) and the 
Gaussian curvature of the surface, and find the leading correction terms in the l/N power 
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series for (Dk(N)) as functions of the curvature. This results in the topological invariance at 
0(1/N). We discuss higher order terms as well, and the case of higher dimensions. Finally, 
we show that a Regge calculus approach provides a simple means of obtaining this topological 
invariance for the case of polyhedral (non-smooth) surfaces. 

2 Preliminaries 

Take any point x, and consider P[Z}fc(iV) = Z] x , the probability density that the point's kth- 
nearest neighboring site lies at a distance / from it. This is equal to the probability density 
of having k — 1 (out of N) sites within distance I, one site (out of N — k + 1) at distance 
/, and the remaining N — k sites beyond distance I. Let us choose units so that our surface 
has total area 1. Since sites are distributed uniformly over the surface, the probability of a 
site lying within distance / is then simply the area A(l) K of a disc of radius / about point x 
on the surface. This is shown in Figure |1| Dropping the argument x (in order to simplify 
the notation), we may then write 




N-k 



Figure 1: A point x and its neighborhood within a distance /. Area A(l) x is the probability 
of finding a site within distance I of x. Here, k — 1 sites lie within this region, the /cth-nearest 
site lies at distance I exactly, and the remaining N — k sites lie outside the region. 
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giving the expectation value (first moment) 

POO 

(D k (N)} = / P[D k {N)=l\ldl 



ri[A(l)]^[l-A(l)r k ^ldL 



(N-k)\(k-l)\ Jo WJ dl 

Using the variable transformation w = A(l), this may be written in terms of the inverse 
function A~ 1 (w) as 

/V! r 1 

ly ■ / A-lf...\ ...k-l/i „..\N—k 



W» = pvlWliji /. A W 10 (1 - "> 

If A (it?) admits the power series expansion in w. 



A = u; 7 ^ Cjto- 7 , for some 7 G [0, 1), (1) 

j=0 

then 

wa, » = Ji* -in g'»/;«^- , (i-») w -*«»»- <2) 

Recognizing the integral as the Beta function B(fc+j+7, N — k + 1) = (k+j+j— 1)! (N — k)\/ 
(iV + j+7)!, 

/no\m (fc + j+7-1)! iV! r . 

Several comments are in order concerning (Dj.(N)). First of all, although we restrict 
ourselves to discussing the first moment of Dk(N), we could in fact consider any moment 
(D k (N)) by taking [A _1 (w)] a instead of A~ l (w) in ([!]). Doing so would alter 7 and the c/s, 
but would not change our results qualitatively. Second of all, there is no loss of generality in 
taking our total surface area to be unity; scaling this area by a constant (or even, as might 
be more intuitive in statistical physics, by N) would provide only a trivial scaling factor in 
our results. Third of all, we could imagine that the point x we consider is itself an (iV + l)th 
site. This is simply a question of nomenclature: the problem of finding the expected distance 
from an arbitrary point to its kth-nearest site, for a system of iV sites, is equivalent to the 
problem of finding the expected distance between kth-nearest neighboring sites, for a system 
of iV + 1 sites. 

We now turn to the properties of A~ 1 (w), and their consequences on (Dk(N)). 
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3 Flat Surfaces 



On a flat surface, if we could neglect edge effects, the area included within distance I would 
simply be A(l) = nl 2 . In that case, we would have A~ l {w) = ^Jw/n, and so from ([!]) and 

(D, 

(D (N)) - 1 itlll^l Nl (A) 

{Dk[N)) -^ (fc-i)! (ivTW ( ) 



There would thus be a complete separation of the /c-dependence and the iV-dependence. 

As we are working with a surface of fixed (unit) area, however, we cannot avoid con- 
sidering edge or finite size effects. Let us restrict ourselves to the case where the surface is 
everywhere locally Euclidean within some minimum neighborhood of radius Iq > 0. (The sim- 
plest example of this is a unit square with periodic boundary conditions, for which l = 1/2. 
Clearly, many other constructions are possible.) Any required modification to the A^ 1 {w) 
expression in flU) then concerns only w greater than w = A(l ). Correspondingly, (§) re- 
mains valid, up to remainder terms from the region of integration wq < w < 1. Since the 
(1 — w) N ~ k ~ 1 term in the integral is bounded above by (1 — Wq) n ~ 1 within this region, 
these remainder terms are exponentially small in N. Equation is thus still correct to all 
orders in a 1/N series expansion, and may be written as 



wm- 1 ( *~ 1/2)! 1 



3 / i \ 

(Jfe-1)! s/N L 1 _ 8iV + ° \PP ) 

where all orders in the series are independent of k. We therefore see that the large iV scaling 
law for /cth-nearest neighbor distances on a 2-D flat surface without a boundary exhibits the 
universality in k to all orders in 1/N. 

The same holds true for flat manifolds in any dimension d. We assume there is some lo 
such that the volume included within distance I < lo is simply the volume of a (^-dimensional 
ball: 

id 



1 

7T 



W 



1/d 



A-\w) = 

As before, this condition allow us to write (D^(N)) up to remainder terms that are expo- 
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(a) I represents chord length distance 



(b) I represents arc length distance 



Figure 2: A spherical surface, where A(l) denotes the area of a spherical cap spanned by /. 
In (a), I represents the chord length distance; in (b), I represents the arc length distance. 

nentially small in N, so from (H), 



Thus for flat spaces without a boundary, of any dimension d, the universality in k holds to 
all orders in 1/N: the fc-dependence and the iV-dependence separate. 

It may be interesting to consider a slight variation on the problem, giving this universality 
exactly and not only to all orders. Take the case of a spherical surface embedded in 3-D 
Euclidean space, with the usual measure of area over the sphere, but with a peculiar sort 
of "distance": rather than the conventional choice of the arc length (geodesic) metric, use 
the chord length. (See Figure |2]a.) For a chord of length I originating at a pole of the 
sphere, the area of the spherical cap spanned by it is simply A(l) = irl 2 . The /cth-nearest 
neighbor distance properties using chord length "distance" on this curved surface then appear 
analogous to those on a flat surface. There is, however, one important distinction. The 
relevant threshold wo for edge effects is in this case wo = 7r(2i?) 2 , where R is the radius of 
the sphere. Since ir(2R) 2 is exactly equal to the total surface area of the sphere, it is set to 
1. Equation (fj) thus requires no corrections at all, and so the universality in @ is exact. 



( D k{N)) 



[(d/2)\} 1/d (Jfe-l + l/d)! N\ 



v^F (k — 1)! (N + l/d)\ 



P/2)!r (fc-l + l/d)! N - 1/d [ 
A (k-l)\ [ 



1/d+l/d 2 
2N 
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4 Curved Surfaces 



Now consider the case of a surface with intrinsic curvature, with the distance defined in 
terms of a metric, i.e., along geodesies of the surface. Let us begin with a spherical surface, 
this time letting I represent arc length; the area of the spherical cap spanned by an arc 
originating at a pole of the sphere (see Figure 0b) is given by 

I 



A(l) 



sphere 



2nR z 



1 — cos ■ 



R 



AtiR 2 sin 2 



I 

2R' 



If the total surface area, 4irR 2 , is normalized to 1 

ere 

A (w) 

sphere 



sin 



so 



sin 



w 



71 

oo 



(5) 



As in the case of the chord length "distance", this A 1 (w) expression is exact everywhere 
for < w < 1. Equations (|l|) and @ then require no corrections, and we find 

1 ^ 1 (2j)\ (k + j-l/2)\ y/NN\ 



(D k (N)) sphe 



^ § 2j + 1 2V (j!)2 (k - 1)! (N + j + 1/2)! 



3=0 

1 (k — 1/2)! 

7T (k 



1 + 



Ak-7 
24N 



(6) 



1)! VN 

Clearly, the k universality does not apply here: the 0(1 /N) coefficient explicitly contains k. 

Another sort of universality, however, is found when we turn to the more general case 
of an arbitrary closed surface, i.e., an abstract 2-D manifold with no boundary. Given a 
smooth surface, we may introduce a system of curvilinear coordinates u and v and write (at 
least piecewise on the surface) the differential length element ds in the conformal, orthogonal 
form [|J: 

ds 2 = f(u,v) [du 2 + dv 2 }. (7) 

The Gaussian curvature K(u, v) of the surface is then expressed in terms of the function 
f(u,v) by 



K 



1 

2/3 



'df 



df 



,d 2 f £ d 2 f 



du 



dv 



du 2 



f 



dv" 1 
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What is A(l) on this surface? To find out, we first determine the manifold's geodesic lines. 
For ds given by (0), we may use the geodesic equation ||: 



d 2 u J_ df_ 



' du 
ds 



' dv V 
ds 



1 df du dv 
H = 0. 

/ dv ds ds 



(9) 



Let us expand u and v as functions of distance s from an initial point, along a fixed geodesic: 

s 2 

u(s) = un + su'n + —u'n + ■ ■ ■ and 



v(s) 



Wo + su'q + —Uq + ■ ■ 
s 2 

v o + sv' + —v'o + ■ ■ 



(10) 



where Uq = u(0), u' = u'(0), etc., and likewise for v. Then, expanding f(u,v) in terms of u 
and v and substituting (|ToD , 



f(u,v) = f(u ,v ) + s[u' f u (u ,v ) + v' f u (u ,v )} + s 2 

(u') 2 (v 1 ) 2 



u" v" 

Y fu{UQ, V ) + y f v (u , V ) 



+ 0(s 3 ), 



where subscripts on / denote partial derivatives. 

Using (0) and (H), we can solve for all but three of the coefficients in ([U]). Let us choose 
Uq, Vo and u' Q to be these three. Now, consider the area A{1) about the point (u ,v ). For 
ds given in (|7|), the differential surface element will be dfi = f dudv, so: 



A{1) 



f du dv 

U, V 



If 



du dv 



ds du'n 



du dv 



ds duL 



11) 



ds du'o du'o ds 

The limits of integration over s are and /; the limits of integration over u' Q , which may be 



found from (0), are —yl/f(uo,Vo) and yl/f(uo,Vo). In evaluating the Jacobian some care 
must be taken, as a sign ambiguity allows two solutions for the coefficients in (PH[). A(l) will 
be the sum of ([□]) evaluated at each of the two solutions, ultimately causing all odd powers 
of I to vanish. 

The result, after lengthy algebraic manipulations, may be written as 



A(l) = ttV 



i - ^ js (fu + il - - fU + o(i 4 



(12) 
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where in order to avoid cluttering the notation, we have omitted the (u ,v ) arguments at 
which all functions are to be evaluated. For the leading correction term in A(l) given in 
(|T2|), we recognize the expression (|8p for the Gaussian curvature K. In retrospect, this is 
not surprising: by symmetry, the correction series to ttI 2 can contain only even powers of 
and if we consider A{1) as a geometric expansion about a flat space approximation, K will 
be the only scalar curvature quantity with dimensions of l~ 2 0. 

With some perseverance, one may carry the expansion in (|T2|) to higher orders, obtaining 



A(l) 



^ _ i_ j \ fnr.-2 or-' 

12 720 



[2K 2 - 3V 2 K) 



I 6 



[8K 3 - 3[10(ViT) 2 + UKV 2 K - 5V A K}) + 0(/ 8 ) 



(13) 



161280 

where V is the gradient operator. We thus obtain a series expansion giving the area of a 
disc on a smooth 2-D surface, in a form that depends only on intrinsic quantities, i.e., not 
on the choice of coordinate system. 

We may invert ( |T3"D to obtain the power series 



w 



K 9K 2 + AV 2 K 
1 + + 



24vr 



1920vr 2 



vr H „ w s 



21504vr 3 



(14) 



As an example, take the special case of a spherical surface, where the Gaussian curvature 
is a constant K = 1/R 2 , or K = Att for a unit surface. All derivatives of K then vanish, 
leaving 

A-\w) s 



/sphere 



w 3w 2 5w 3 4 , 

H 1 1 h 0{w A ) 

6 40 112 v ; 



from which we recover the first few terms of our earlier result 

Given an expression for A~ 1 (w) on a general 2-D surface, we may now find (Dj.(N)) using 
(jl|) and (0). It is helpful at this point to define the reduced variable (Dk(N)) by dividing 
out the leading asymptotic (large N) behavior from (Dk(N)). Recalling the notation of (|T|), 
A^ 1 (w) = w 1 J2j c j w ^ f° r some 7 e [0, 1), define: 

(D k (N)} = (D k (N)} — 



c (fc + 7- 1)! 



AT7 



Arr y C J (fc + J+7-1)! 

^ c (fc + 7-1)! (iV + J+ 7)!' 



(15) 



j=0 



so that lim7v^oo(-Dfc(^)) = 1 (this is a consequence of Stirling's law). (D k (N)) then corre- 
sponds to the 1/N power series that we have frequently seen appearing in (Dk(N)), giving 
the corrections to its leading (large N) behavior. It is in this power series that the inter- 
esting universality properties emerge. Consider the average / (D^(N)) dfi over the entire 
surface, obtained — via ([[]) and (|l^) — from precisely the average of the series coefficients 
in ([14]) . Examine, in particular, the 0(w) term of ([14]). By the Gauss-Bonnet theorem [[|, 
/ K dfi = 2ttx on any closed surface, where \ is the Euler characteristic of the surface, a 
topological invariant. Up to leading corrections, then, 



A 1 (w)d/2 
(D k {N))dfi 



1 + -^w + 0(w 2 ) 
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giving 



NN\ 



(N + 1/2)! 



12iV + 3/2 \N 2 



x (2fc+l)-9 1 
X + 7777 +° m ' ( 16 ) 



24iV ViV 2 



We thus discover a different sort of universality from the one we had in the case of flat space. 
To 0(1/N), the scaling law for kth- nearest neighbor distances depends only on the surface's 
topology, and not on its detailed properties. 

The Euler characteristic x for a surface is related to its genus g by x — 2(1— g) . Taking the 
torus as one example, g = 1, so \ = an d the /c-dependence in (^) once again disappears, 
at least at 0(1/N). This is to be expected: a flat space with periodic boundary conditions 
has, after all, the topology of a torus. And conversely, because of the topological invariance, 
all tori behave like flat space to 0(1/N). Taking the spherical surface as another example, 
g = 0, so x = 2 and we recover from (|16D the power series in @. 

The properties of J{D/,(N)) dfi are far less clear at higher orders in 1/N. Using the 
divergence theorem and integration by parts, we may obtain from (|14]) 



A 1 (w)dfi 



K 3K 2 2 15K 3 + 16KV 2 K 3 _ 4N 

1 H to H =■ w H ^ wr + 0(w 4 

24tt 640tt 2 21504tt 3 v 7 



dfi. (17) 



If we looked only at terms up through 0(w 2 ), we might believe that this series is simply, by 



analogy with (j|), the expansion of f(2/y/K)sin~ 1 yj Kw/An dfi. Unfortunately, starting at 
0(w 3 ) we see this is not true, since the contributions of curvature and its gradients do not 
all vanish in the average over the surface! Furthermore, even for terms in ([17]) of the form 
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/ K n dfi, at n > 1 there is no straightforward equivalent to the Gauss-Bonnet theorem; the 
theorem is a direct consequence of the integrand's linearity. Thus for a general 2-D surface, 
a simplified form does not appear to exist for the terms in J{Dk(N))dfj, beyond 0(1/N). 
More particularly, the only case in which / (Dk(N)) d\i would be independent of k beyond 
0(1 /N) is if the curvature is identically equal to 0, i.e., a flat surface. 

Let us briefly consider the case of curved higher-dimensional manifolds. The calculation is 
now far more complicated, as it is no longer possible to write the metric tensor in a conformal 
form as we did in fl7|). In addition, whereas in 2-D the only intrinsic scalar quantity describing 
curvature is the Gaussian curvature K, for d > 2 there are d(d — — 2)(d + 3)/12 different 
such quantities [Q . However, all of them except K itself have dimensions of higher order than 
l~ 2 . It thus seems reasonable to conjecture that, as we argued in 2-D, the 0(l 2 ) correction 
term in A(l) can only involve K. (Indeed, we have verified that this is true in 3-D.) In 
that case, we may rely on the example of the spherical surface — easily generalized to d 
dimensions — to provide us with the initial terms for a general manifold: 

^12 



A(l) 
A-\w) 



r 



(dp) 
(d/2)\ 



d(d-l) K 
6~ 



d 



or 



-i l/d 



7T 



d/2 



W 



l/d 



1 + 



d—\K /(d/2)!' 
d + 2~Q I 7r d / 2 



2/d 



w 2/d + 0{w A/d ) 



Note that A 1 (w) now contains a series in w 2 l d rather than in w. Appropriately modifying 
([]]), it may then be shown that (Dk(N)) is in general given by a series in for odd d, 

and l/N 2 > d for even d. 

Consider, finally, the average / A (w) dfi over the manifold. The higher-dimensional 
generalization of the Gauss-Bonnet theorem || involves an integrand of 0(l/l d ), or 0(l/w). 
The leading correction term from (|I8|), / K dfi, therefore cannot be simplified further for 
d > 2; the only term that could possibly give rise to a topological invariant is the coefficient 
at 0(l d ), or 0(w). If d is odd, it is rather certain that no topological invariant will exist in the 
series. If d is even, the 0(w) term will first contribute to the / (Dk(N)) d\x series at 0(1/N) 
- as in 2-D, although at higher dimensions this will no longer be the leading correction term. 
While one cannot rule out the possibility of obtaining a topological invariant at 0(1/N), 
the 0(w) term in A~ 1 (w) is in general a complicated one involving many different curvature 
scalars, and so this is far from obvious. We leave it as an open question. 
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5 Regge Calculus 



We have remarked that from a physical point of view it is natural, in the 2-D case, for 
the leading corrections in A(l) to contain only the Gaussian curvature, as this represents 
the leading deviation from planarity. Consequently, only the mean curvature — or, using 
the Gauss-Bonnet theorem, the Euler characteristic \ — matters in the 0(1/N) term of 
/ (Dk(N)) dfi. We have seen using differential methods (geodesies) that this physical picture 
is indeed correct. These methods apply to a smooth closed surface. For polyhedral surfaces, 
which are not smooth, we may in fact obtain a similar result more easily, using the non- 
differential method of Regge calculus. Consider a polyhedron with a number of vertices, 
edges, and faces. Following the work of Regge and others since then |7[ , we observe that 
the curvature is concentrated at the vertices and is measured by a deficit angle: if 9i is the 
sum of the angles incident on vertex i, the deficit angle at that vertex is A, = 27r — 0{. It may 
then be shown that the Gauss-Bonnet theorem, on polyhedra, reduces to Euler's relation 
27cx = E i A i . 

Let P be a polyhedron with a fixed number of vertices, and consider the problem of 
finding the large N scaling series / (Dk(N)) dfi on P. As N — > oo, corrections to the flat 
space value about a given point x arise only when x is near one of the vertices, because 
only in that case can curvature (i.e., the deficit angle) enter into the local calculation of A(l) 
about x. It is then sufficient to understand the corrections associated with one vertex at a 
time. Consider the neighborhood of a vertex i. A(l) receives a correction from the flat space 
value, and by a simple geometric construction, one can see that this correction is exactly 
proportional to the deficit angle Aj. Correspondingly, the leading correction term both in 
A^ 1 (w) and in the scaling series will be proportional to Aj for small deficit angles. Now sum 
over all the vertices % — 1, . . . , N, assuming that all the deficit angles are small. We then find 
that the 0(1 /N) term in J (Dk(N)) d[i is proportional to %i an d we recover the topological 
invariant derived in the case of a smooth manifold. 

A word of caution is necessary, however. It is tempting at this point to take the limit 
where P becomes a smooth manifold, expecting to recover flTEp. Unfortunately this will not 
work; a direct computation shows that the limit does not commute with the limit N — ► oo 
taken above, and the coefficient thus obtained at 0(1/ N) will not be the correct one. 
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6 Conclusions 

Given iV sites distributed randomly and uniformly on a surface with no boundaries, we have 
considered the properties of mean distances to neighboring sites. When the surface is flat, we 
have seen that in the expression for the mean fcth-nearest site distance, the fc-dependence and 
iV-dependence separate. The scaling law in N for mean /cth-nearest neighbor distances is thus 
independent of k. This universality applies equally well to higher moments of the distances, 
and to Euclidean manifolds in dimensions greater than 2. For surfaces with curvature, while 
this general property is no longer valid, we have found that when the fcth-nearest neighbor 
distance is written as a large N expansion, averaged over the surface, the leading correction 
term in the series is a topological invariant. The scaling series thus depends, to 0(1/N), on 
the genus of the manifold but not on its other properties. 

Although we have considered these universalities only for the moments of point-to-point 
distances, similar properties hold for higher order simplices such clS cirCcLS of triangles as- 
sociated with nearby points. The problem is thus a natural one to consider further in the 
context of random triangulations, foams or other physical problems |2j tightly connected to 
geometry. 
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